Computational Fluid Dynamics Model Research Articles

Two-Fluid and Discrete Population Balance Model for Simulation of Batch Sedimentation Relevant to Plutonium Reconversion

There are widespread occurrence and application of solid-liquid sedimentation processes among different industries. Therefore, it becomes important to understand the hydrodynamic inside the process equipment and particle agglomeration characteristics. In the present work, solid-liquid sedimentation is analyzed, which will be helpful for the design of continuous process equipment in plutonium (Pu) reconversion. Experiments were carried out in a batch settler to understand solid sedimentation in suspension in terms of varying the overall solid fractions. Euler-Euler two-fluid simulations were performed to investigate the local and overall solid phase volume fraction distributions, position of the active interface (AI) (settling curve), axial solid phase velocity, and pressure distributions during settling, and selected data were compared with the measured data. Further, the discrete population balance model (PBM) with different agglomeration kernels was used with the two-fluid computational fluid dynamics (CFD) model to understand and further improve predictions in terms of the AI position. The variation in number density of the different particles in the settler with time was investigated. The predicted results show that agglomeration is dominant during the sedimentation process and application of the discrete PBM with the CFD model enhances the predictive capability in comparison with the predictions obtained from only the two-fluid model. The results reported using CFD+PBM will aid in the design of continuous process equipment (thickener/precipitator) for Pu reconversion.

Insights into reaction–diffusion behaviors of acetylene selective hydrogenation on pellet catalysts

AbstractIn this work, a particle‐resolved computational fluid dynamics model of the acetylene hydrogenation process is developed to investigate the effects of catalyst particle structures on the reaction–diffusion behaviors aiming to improve selectivity toward target ethylene. The effects of packing structures on mass and heat transfer are explored by employing particles with varying shapes, wall thicknesses, and external diameters. The simulation results reveal that decreasing diffusion paths and elevating reactor bed temperature will enhance ethylene selectivity, and cylinder and Raschig ring packing structures exhibit the lowest and highest ethylene selectivity of 35.6% and 48.9% at 70% of acetylene conversion, respectively. Reducing wall thicknesses of Raschig ring particles facilitates the diffusion of generated ethylene from the interior zone of catalysts but concurrently inhibits the conversion of acetylene to ethylene. The Raschig ring catalyst particle with 1.9 mm of wall thickness and 3.5 mm of external diameter is finally revealed to exhibit the highest ethylene selectivity.

Analysis of Anthropogenic Waste Heat Emission from an Academic Data Center

The rapid growth in computing and data transmission has significant energy and environmental implications. While there is considerable interest in waste heat emission and reuse in commercial data centers, opportunities in academic data centers remain largely unexplored. In this study, real-time onsite waste heat data were collected from a typical academic data center and an analysis framework was developed to determine the quality and quantity of waste heat that can be contained for reuse. In the absence of a comprehensive computer room monitoring system, real-time thermal data were collected from the data center using two arrays of thermometers and thermo-anemometers in the server room. Additionally, a computational fluid dynamics model was used to simulate temperature distribution and identify “hot spots” in the server room. By simulating modification of the server room with a hot air containment system, the return air temperature increased from 23 to 46 °C and the annual waste heat energy increased from 377 to 2004 MWh. Our study emphasizes the importance of containing waste heat so that it can be available for reuse, and also, that reusing the waste heat has value in not releasing it to the environment.

Numerical Analysis of Air Supply Alternatives for Forced-Air Precooling of Agricultural Produce

Precooling agricultural produce is an intensive, energy-consuming process. To improve the efficiency of forced-air precooling and ultimately contribute to energy sustainability for postharvest storage of fresh produce, we designed three alternative air supply systems, simulated their cooling performances over a 96 h precooling process in a cold storage facility storing Chinese cabbages, and then compared their performances with a conventional design. All models were developed on a large scale on the basis of validated computational fluid dynamics models. The horizontal air supply scheme shortened the seven-eighths cooling time by 18.8%, and its maximum cooling rate increased by 19.7% compared to the conventional air supply scheme. The seven-eighths cooling time under another alternative design, the vertical air supply scheme, was 9.4% lower than the conventional, with the maximum cooling rate increasing by 10.5%. However, the maximum cooling rate of the last alternative design, the perforated ceiling air supply system, was 6.6% less than the conventional scheme, resulting in a 6.3% longer seven-eighths cooling time. The heterogeneity index of temperature implied that the horizontal air supply offered better overall cooling uniformity than the other designs, which can be attributed to its evenly distributed airflows and well-organized air movement paths, based on the combined analysis of temperature contours and air velocity contours at selected planes. Our findings are expected to provide practical guidelines for the refinement of the air supply system to improve its energy sustainability in forced-air precooling.

Evaluating the accuracy of cerebrovascular computational fluid dynamics modeling through time-resolved experimental validation

Accurate modeling of cerebral hemodynamics is crucial for better understanding the hemodynamics of stroke, for which computational fluid dynamics (CFD) modeling is a viable tool to obtain information. However, a comprehensive study on the accuracy of cerebrovascular CFD models including both transient arterial pressures and flows does not exist. This study systematically assessed the accuracy of different outlet boundary conditions (BCs) comparing CFD modeling and an in-vitro experiment. The experimental setup consisted of an anatomical cerebrovascular phantom and high-resolution flow and pressure data acquisition. The CFD model of the same cerebrovascular geometry comprised five sets of stationary and transient BCs including established techniques and a novel BC, the phase modulation approach. The experiment produced physiological hemodynamics consistent with reported clinical results for total cerebral blood flow, inlet pressure, flow distribution, and flow pulsatility indices (PI). The in-silico model instead yielded time-dependent deviations between 19–66% for flows and 6–26% for pressures. For cerebrovascular CFD modeling, it is recommended to avoid stationary outlet pressure BCs, which caused the highest deviations. The Windkessel and the phase modulation BCs provided realistic flow PI values and cerebrovascular pressures, respectively. However, this study shows that the accuracy of current cerebrovascular CFD models is limited.

Simulation of Transient Temperature and Clearance after Shutdown of Aeroengine Based on CFD and FEA Coupled Models

It is crucial to comprehend the heat soak phenomenon, which may result in a significant temperature increase after the shutdown followed by a gradual decrease. This could bring potential risks for the engine including oil coking. The temperature change of engine components dictates the clearance after shutdown, while startup strategies are primarily based on this. A simulation strategy, utilizing computational fluid dynamics (CFD) and finite element analysis (FEA) coupled models, is suggested to investigate the transient temperature and clearance after shutdown. The maximum temperature deviation between the simulation result and experimental data are less than 6%. Flow parameters, including velocity and mass flow rate obtained from the CFD result, were applied as boundaries of the FEA model. Based on the FEA model, transient temperature calculations were also conducted for 20 hours after shutdown. The results indicate that the FEA model demonstrates good agreement with the CFD simulation, with a maximum deviation of less than 5% and at only 0.2% of the simulation time. After the engine shuts down, the stator’s temperature change rate is faster than that of the rotor due to better cooling conditions and relatively small heat capacity. Consequently, the seal clearance increases in the initial period after shutdown and then decreases to a minimum value. The nondimensional minimum clearance can be 0.8 times the cold state value at the location of the high-pressure turbine seal.

Aerodynamic performance and wind-induced response of carbon fiber-reinforced polymer (CFRP) cables

To study the aerodynamic performance and wind-induced response of carbon fiber-reinforced polymer (CFRP) cables, CFRP cable was designed by replacing a steel cable in a tied arch bridge based on stiffness, strength and area equivalent criteria, respectively. The aerodynamic performance and wind-induced response of CFRP cable and steel cable were studied and compared by computational fluid dynamics (CFD) model. Based on the computational results, optimal cable replacement criterion was proposed for CFRP cable to replace steel cable. In addition, surface modification was conducted by engrooving vertical symmetric (VS), vertical asymmetric (VA) and helical symmetric (HS) V-shaped grooves to improve the aerodynamic performance and wind-induced response of CFRP cable. Results showed that CFRP cables exhibited inferior aerodynamic performance and wind-induced response in most cases. However, CFRP cable based on stiffness equivalent criterion exhibited better aerodynamic performance and wind-induced vibration properties compared to the other two cable replacement criteria, thus is regarded as the optimal substitute for steel cable. In addition, HS grooves generated symmetric disturbances and caused approximately equivalent boundary layer separation delays uniformly and continuously along the cable length, thus exhibiting better effect in decreasing the reverse flow region, the maximum negative flow velocity, the vortex shedding frequency and the wind-induced vibration amplitude of CFRP cable. Hence, the stiffness equivalent criterion combined with surface modification with HS V-shaped grooves was proposed to replace steel cable with CFRP cable. This study can provide insights into the aerodynamic performance and wind-induced response of CFRP cable and instructions for cable replacement practice.

A Comprehensive Numerical Study of a Wedge-Shaped Textured Convergent Oil Film Gap

The modification of surface geometries to reduce friction is an omnipresent topic of research. In nature, different low-friction surfaces, such as fish skins, exist. To transfer this knowledge to technical applications, for example, to journal or plain bearings, many numerical and experimental studies of textured surfaces have been performed. In this work, the influence of the geometric parameters (texture length l, width b, angle α and start position xstart) of a wedge-shaped texture on three different convergent oil film gaps was analyzed in full-film lubrication and compared with untextured oil film gaps. With the aid of a CFD (computational fluid dynamics) model, a comprehensive variation study was conducted, and the best-performing wedge-shaped texture was determined. The results show that an open texture at the inlet provides the largest improvement. Furthermore, it can be observed that the optimal relative texture width and absolute inlet height for the three investigated oil film gaps are similar. In contrast to the volume flow of the untextured geometry, the volume flow of the textured one is significantly higher, especially that perpendicular to the movement direction.

Influence of building spatial patterns on wind environment and air pollution dispersion inside an industrial park based on CFD simulation.

The "spatial pattern-wind environment-air pollution" within building clusters is closely interconnected, where different spatial pattern parameters may have varying degrees of impact on the wind environment and pollutant dispersion. Due to the complex spatial structure within industrial parks, this complexity may lead to the accumulation and retention of air pollutants within the parks. Therefore, to alleviate the air pollution situation in industrial parks in China and achieve the circular transformation and construction of parks, this study takes Hefei Circular Economy Demonstration Park as the research object. The microscale Fluent model in computational fluid dynamics (CFD) is used to finely simulate the wind flow field and the diffusion process of pollutants within the park. The study analyzes the triad relationship and influence mechanism of "spatial pattern-wind environment-air pollution" within the park and studies the influence of different spatial pattern parameters on the migration and diffusion of pollutants. The results show a significant negative correlation between the content of pollutants and wind speed inside the industrial park. The better the wind conditions, the higher the air quality. The spatial morphology parameters of the building complex are the main influences on the condition of its internal wind environment. Building coverage ratio and degree of enclosure have a significant negative correlation with wind conditions. Maintaining them near 0.23 and 0.37, respectively, is favorable to the quality of the surrounding environment. Moreover, the average height of the building is positively correlated with the wind environment condition. The rate of transport and dissipation of pollutants gradually increases as the average building height reaches 16 m. Therefore, a reasonable building planning strategy and arrangement layout can effectively improve the wind environment condition inside the park, thus alleviating the pollutant retention situation. The obtained results serve as a theoretical foundation for optimizing morphological structure design within urban industrial parks.

A critical analysis of the influence of architecture on the temperature field of RC structures subjected to fire using CFD and FEA models

The temperature–time-heating curve describes the dynamics of a fire in a room. A variety of factors related to the architecture and the thermal properties of the furniture influence the temperatures. The reinforced concrete (RC) structures of the buildings are designed for the case of fire with standardized curves that ignore the real fire temperatures. In this study, a group of individual rooms were solved numerically to understand the thermal cross-sectional field in RC structures based on real fire curves. The fire dynamics of the room were defined by a computational fluid dynamics (CFD) model solved with the Fire Dynamics Simulator (FDS) software. A finite element analysis (FEA) model was created for the RC structures and solved using Abaqus software. A typical RC ribbed slab on the roof of the room was considered. The structure was heated according to the real and standardized fire curves. The real fire temperatures were up to 70 % lower than those predicted by the standard. The results show how safe the buildings are from fire-induced collapse.

GPU-accelerated CFD-DEM modeling of gas-solid flow with complex geometry and an application to raceway dynamics in industry-scale blast furnaces

The coupling of Computational Fluid Dynamics (CFD) and Discrete Element Model (DEM) is a powerful tool for simulating dense particulate systems, yet the conventional CFD-DEM has limits for systems with large particle numbers and complex geometry. This paper reports a novel GPU-based CFD-DEM model to simulate the gas–solid flow with large particle numbers and complex geometry. A novel coupling strategy between the CFD solver and DEM solver is developed, featuring high efficiency and stability. The developed model is validated against the experimental measurements, and its efficiency is compared to the previous CFD-DEM simulations. Then, for demonstration, the model is employed to simulate the dynamic behavior of gas–solid flow in the raceway in ironmaking blast furnaces by considering complex tuyere structure details and the huge particle numbers involved. This model allows to study the effect of the tuyere angle in terms of raceway formulation and tuyere erosion. The results show that the largest and most stable raceway volume can be reached at 5° downward tuyere, although the −5° tuyere nose experiences more wear than 10° downward tuyere. The model provides a cost-effective tool to overcome the longstanding challenge of simulating dense fluid-particle systems with huge particle numbers and complex geometry.

Impact of Age-Related Change in Caval Flow Ratio on Hepatic Flow Distribution in the Fontan Circulation.

The Fontan operation is a palliative technique for patients born with single ventricle heart disease. The superior vena cava (SVC), inferior vena cava (IVC), and hepatic veins are connected to the pulmonary arteries in a total cavopulmonary connection by an extracardiac conduit or a lateral tunnel connection. A balanced hepatic flow distribution (HFD) to both lungs is essential to prevent pulmonary arteriovenous malformations and cyanosis. HFD is highly dependent on the local hemodynamics. The effect of age-related changes in caval inflows on HFD was evaluated using cardiac magnetic resonance data and patient-specific computational fluid dynamics modeling. SVC and IVC flow from 414 patients with Fontan were collected to establish a relationship between SVC:IVC flow ratio and age. Computational fluid dynamics modeling was performed in 60 (30 extracardiac and 30 lateral tunnel) patient models to quantify the HFD that corresponded to patient ages of 3, 8, and 15 years, respectively. SVC:IVC flow ratio inverted at ≈8 years of age, indicating a clear shift to lower body flow predominance. Our data showed that variation of HFD in response to age-related changes in caval inflows (SVC:IVC, 2, 1, and 0.5 corresponded to ages, 3, 8, and 15+, respectively) was not significant for extracardiac but statistically significant for lateral tunnel cohorts. For all 3 caval inflow ratios, a positive correlation existed between the IVC flow distribution to both the lungs and the HFD. However, as the SVC:IVC ratio changed from 2 to 0.5 (age, 3-15+) years, the correlation's strength decreased from 0.87 to 0.64, due to potential flow perturbation as IVC flow momentum increased. Our analysis provided quantitative insights into the impact of the changing caval inflows on Fontan's long-term HFD, highlighting the importance of SVC:IVC variations over time on Fontan's long-term hemodynamics. These findings broaden our understanding of Fontan hemodynamics and patient outcomes.

INFLUENCE OF SPRAY DISTANCE ON THE POROSITY OF Ni-BASED AMORPHOUS COATINGS: NUMERICAL SIMULATION AND EXPERIMENT

The influence of injection distance on the porosity of Ni-based amorphous coatings (AMCs) prepared with a high-velocity air fuel (HVAF) process is discussed based on a numerical analysis and experimental methods. A computational fluid dynamics model was established to demonstrate the gas flow field and behavior of particles in flight at different spraying distances during HVAF spraying. When analyzing the changes in the particle velocity and temperature, the spraying distance is less than 30 µm. The velocity and temperature changes of small particles have a significant impact, and the optimal spray distance (350 mm) for obtaining a low porosity coating is predicted. The calculation was validated experimentally by producing a Ni-based AMC with a low porosity (1.87 %) that was manufactured using the predicted HVAF optimal spraying distance.

Machine learning enhanced exploration of bubble dynamics beneath a horizontal wall

A better understanding of the dynamic behavior of a single bubble rising and impacting on a horizontal substrate is crucial, yet it has traditionally been constrained by the necessity for demanding experiments and extensive time-consuming simulations. This work develops a machine learning-enhanced model to determine characteristic quantities for the dynamic behavior of rising bubbles beneath the horizontal wall based on Computational Fluid Dynamics (CFD) simulation. We establish an unprecedentedly extensive dataset from accurate CFD modeling, covering a vast range of values for input parameters: Galilei number (1∼100), Bond number (0.001∼20), and rising distance (0.1R∼40R). Three popular regression models (Support Vector Regression, Back Propagation Neural Network, and Random Forest) are trained to capture key features for the first bounce of a bubble after collision, exhibiting promising predictive capacity, generalization, and interpretability. Our current approach facilitates the generation of dense datasets, enabling finer interpolation and detection of the physical properties of the unknown fluids. This pioneering data-driven methodology paves the way for advancing our understanding and prediction of bubble dynamics, marking a significant contribution to fluid dynamics research.

3D concrete printing using computational fluid dynamics: Modeling of material extrusion with slip boundaries

This paper investigates the role of slip boundary conditions in computational fluid dynamics modeling of material extrusion and layer deposition during 3D concrete printing. The mortar flow governed by the Navier-Stokes equations was simulated for two different slip boundary conditions at the extrusion nozzle wall: no-slip and free-slip. The simulations were conducted with two constitutive models: a generalized Newtonian fluid model and an elasto-viscoplastic fluid model. The cross-sectional shapes of up to three printed layers were compared to the experimental results from literature for different geometrical- and speed-ratios. The results reveal that employing free-slip boundary conditions at the extrusion nozzle wall improves layer-mimicking quality for both constitutive models, indicating the presence and importance of a lubricating layer of fine particles at the concrete-solid wall interface. This enhanced performance is primarily due to the observed decrease in extrusion pressure that minimizes layer height- and width-deviations compared to the experimental prints. Furthermore, the free-slip boundary conditions play an important role in predicting the multilayer prints, its deformation and groove shapes.

Particle separation mechanisms in suspension-feeding fishes: key questions and future directions

Key unresolved questions about particle separation mechanisms in suspension-feeding fishes are identified and discussed, focusing on areas with the potential for substantial future discovery. The published hypotheses that are explored have broad applicability to biological filtration and bioinspired improvements in commercial and industrial crossflow microfiltration processes and microfluidics. As the first synthesis of the primary literature on the particle separation mechanisms of marine, estuarine, and freshwater suspension-feeding fishes, the goals are to enable comparisons with invertebrate suspension-feeding processes, stimulate future theoretical and empirical studies, and further the development of biomimetic physical and computational fluid dynamics models. Of the eight particle separation mechanisms in suspension-feeding fishes, six have been proposed within the past twenty years (inertial lift and shear-induced migration, reduction of effective gap size by vortices, cross-step filtration, vortical flow along outer faces of gill raker plates, ricochet filtration, and lateral displacement). The pace of discovery is anticipated to continue accelerating. Multidisciplinary collaboration and integration among biologists and engineers (including chemical, mechanical, biomedical, and filtration engineering) will result in new perspectives to identify patterns and potential unifying mechanisms across the breadth of suspension-feeding fish taxa, morphology, and function.

Modeling electrolysis in reduced gravity: producing oxygen from in-situ resources at the moon and beyond

Molten Regolith Electrolysis, as an in situ resource utilization (ISRU) technology, has the potential to enable the production of oxygen and metallic alloys on the Lunar surface; opening new doors in Cis-Lunar, and eventually Martian space exploration. This research studies the fundamental physics which govern the formation, growth, detachment, and rise of electrolytic bubbles. To this end, computational fluid dynamic (CFD) models were developed and run, to simulate water electrolysis, molten salt electrolysis (MSE), and molten Lunar regolith (MRE) electrolysis across multiple reduced gravity levels. The results demonstrate that reduced gravity, electrode surface roughness (possibly due to surface degradation), fluid properties, and electrode orientation can all affect electrolytic efficiency and possibly even stall electrolysis by delaying bubble detachment. The findings of this research must be considered when designing and operating electrolysis systems at reduced gravity levels.

A comprehensive model for viscoplastic flows in channels with a patterned wall: longitudinal, transverse and oblique flows

We develop a comprehensive model for the creeping Poiseuille Bingham flow in channels equipped with a patterned wall, i.e. decorated with grooves or stripes that may represent a superhydrophobic (SH) or a chemically patterned (CP) surface, respectively, with longitudinal, transverse and oblique groove (stripe) orientations with respect to the applied pressure gradient. We rely on the Navier slip law to model the boundary condition on the slippery grooves. We develop semi-analytical, explicit-form and complementary computational fluid dynamics models, with solutions that have reasonable agreement. In contrast to its Newtonian analogue, a distinct solution for the oblique configuration, with an a priori unknown transform matrix, must be developed due to the viscoplastic nonlinear rheology. Our focus is to systematically analyse the effects of the Bingham number ( $B$ ), slip number ( $b$ ), groove periodicity length ( $\ell$ ), slip area fraction ( $\varphi$ ) and groove orientation angle ( $\theta$ ), on the slip velocities, effective slip length ( $\chi$ ), slip angle difference ( $\theta -s$ ), mixing index ( $I_M$ ), flow anisotropy and flow regimes. In particular, we demonstrate that, as $B$ increases, the maximum values of the shear component of $\chi$ , $\theta -s$ and $I_M$ occur progressively at smaller values of $\theta$ , compared with their Newtonian counterparts.

Simulation of liquid cone formation on the tip apex of indium field emission electric propulsion thrusters

Field emission electric propulsion (FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities, making them ideal propulsion devices for micro-nano satellites. However, the detection of certain aspects, such as the evolution process of the liquid cone and the physical quantities at the cone apex, proves challenging due to the minute size of the needle tip and the vacuum environment in which they operate. Consequently, this paper introduces a computational fluid dynamics (CFD) model to gain insight into the formation process of the liquid cone on the tip apex of indium FEEP. The CFD model is based on electrohydrodynamic (EHD) equations and the volume of fluid (VOF) method. The entire cone formation process can be divided into three stages, and the time-dependent characteristics of the physical quantities at the cone apex are investigated. The influences of film thickness, apex radius size and applied voltage are compared. The results indicate a gradual increase in the values of electrostatic stress and surface tension stress at the cone apex over an initial period, followed by a rapid escalation within a short duration. Apex configurations featuring a small radius, thick film and high voltage exhibit a propensity for liquid cone formation, and the cone growth time decreases as the film thickness increases. Moreover, some unstable behavior is observed during the cone formation process.

Computational analysis and experiments of spatter transport in a laser powder bed fusion machine

This work focuses on how spatter particles are transported within a laser powder bed fusion (L-PBF) machine. The machine’s gas flow rate and salient flow features are studied with a computational fluid dynamics (CFD) model and are validated with experimental measurements of the flow near the build plate. The CFD model is also paired with a discrete phase model (DPM) to show how a spatter particle’s diameter, speed, ejection angle, material, and ejection location each affect its trajectory. The spatter model results are also validated by experiments wherein hot spatter particles are identified with an infrared camera. Overall, this work brings focus to several spatter-mitigation strategies including increasing gas flow rate, controlling spatter production through laser settings, and changing part placement. Throughout the infrared experiment’s single build cycle, over 14 million spatter particles were detected, contaminating both the parts and the powder bed. Both the spatter model and experiment suggest that spatter particles travel primarily in the downstream direction. This suggests that placing parts next to one another instead of downstream from one another could be a strategy to prevent spatter particles from traveling between parts. A low-speed zone is also discovered beneath the lower nozzle and is predicted to negatively alter the removal of spatter from that region of the build area.